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Abstract 

The main scope of this research consists in evaluating the energy- 
momentum (gravitational field plus matter) and gravitational angular 
momentum densities in the universe with global rotation, considering 
the Godel-Obukhov metric. For this, we use the Hamiltonian for- 
malism of the Teleparallel Equivalent of General Relativity (TEGR), 
which is justified for presenting covariant expressions for the consid- 
ered quantities. We found that the total energy density calculated by 
the TEGR method is in agreement with the results reported by other 
authors in the literature using pseudotensors. The result found for 
the angular momentum density depends on the rotational parameter 
as expected. We also show explicitly the equivalence among the field 
equations of the TEGR and Einstein equations (RG), considering a 
perfect fluid and Godel-Obukhov metric. 
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1 Introduction 



Determination of the energy, momentum, and angular momentum of the 
gravitational field is a long-standing problem in General Relativity (GR). 
In a general way, it is believed that the energy of the gravitational field is 
not localizable, that is, defined in a finite region of space. This problem is 
addressed in the literature by means of different approaches. The most com- 
mon approach is based in the use of pseudotensors, since these gravitational 
physical quantities do not possess proper definitions in terms of tensorial 
equations [1]. 

In order to consider the material content of space-time, it is possible to de- 
fine the energy-momentum (angular momentum) complex as the sum of the 
energy-momentum (angular momentum) pseudotensor of the gravitational 
field and the energy-momentum (angular momentum) tensor of the mat- 
ter. These complexes appear with several names, such as Landau-Lifshitz, 
Bergman-Thompson, and Einstein, among others [2]. They differ from each 
other in the way in which they are constructed. These complexes have been 
applied to several configurations of the gravitational field, as Friedmann- 
Lemaitre-Robertson- Walker (FLRW) [3]-[7], Bianchi I and II [8], [9], Godel 
[10], and Godel-type universes [11] and for a general non-static spherically 
symmetric metric of the Kerr-Schild class [12]. 

In an attempt to contour the localized gravitational energy problem, we 
consider the Weitzenbock space-time, a particular case of Riemann-Cartan 
space-time constrained to have zero curvature. In this space-time, the grav- 
itational field is described in terms of the tetrad field e"^, instead of metric 
tensor gi^j,. In the context of teleparallel gravity, the zero curvature condition 
means that the parallel transport of the orientation of the tetrads is path inde- 
pendent, that is, there is a global moving frame or teleparallelism. Mathemat- 
ically, as a natural consequence of this definition, the Weitzenbock covariant 
derivative of the tetrad field vanishes identically: Vi^^qaj = diyCafi—T'^^eap = 0, 
where F^^ is the Weitzenbock connection given by F^^ = ej^dyt"- y,. The grav- 
itational interaction is attributed to the torsion tensor Ta^u— d^Cau — di^Ca^- 
In Andrade & Pereira [13], the teleparallel gravity can indeed be understood 
as a gauge theory for the translation group. In that approach, the gravita- 
tional interaction is described by a force similar to the Lorentz force equation 
of electrodynamics, with torsion playing the role of force. 

In 1961, M0ller [14] showed that only in terms of tetrads can we obtain a 
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Lagrangian density that leads to a tensor of gravitational energy-momentum. 
This tensor, constructed from the first derivatives of the tetrads, does not 
annul any coordinate's transformation. 

In 2001, Maluf established the Hamiltonian formulation of the Telepar- 
allel Equivalent of General Relativity (TEGR) without the Schwinger time 
gauge. In this formulation, the definition of the energy and gravitational 
angular momentum arises by suitably interpreting the integral forms of the 
constraints' equations C" = and F"^ = 0, respectively [15], [16]. Several 
configurations of gravitational energies were investigated with success, such 
as the space-time configurations of de Sitter [17], conical defects [18], static 
Bondi [19], disclination defects [20], Kerr black hole (irreducible mass) [21], 
BTZ black hole [22], and Kerr anti-de Sitter [23]. The definition of gravi- 
tational angular momentum as the integral form of the constraint equation 
was applied satisfactorily for the gravitational field of a thin, slowly rotating 
mass shell [16] and for the rotational BTZ black hole [24]. 

On the other hand, the idea that the Universe is rotating was first men- 
tioned by Gamow [25]. The exact solution of the Einstein equations for 
the model of a homogeneous, static universe with rotation was proposed 
for Godel [26]. However, the Godel equations presented serious problems, 
such as closed time-like curves (CTCs) that exhibit a causality violation in 
space-time. Several new models based on Godel's original idea were cre- 
ated in an attempt to address the problems in their model [27]-[29]. Among 
these ideas, we stress the nonstationary model proposed by Obukhov [30], 
which presents expansion as well as rotation. This generalization is known as 
the Godel-Obukhov or Godel-type anisotropic cosmological model [31]. This 
model would not be in confiict with any cosmological observations [32]. 

Despite advances in the theoretical understanding of these problems, the 
global rotation has not yet been detected. However, we should stress that 
there are weak observational evidences of global rotation of the universe [33] . 

In this work, we show explicitly the equivalence between GR and TEGR. 
More specifically, we consider the solution for an anisotropic and homoge- 
neous universe described by the Godel-Obukhov metric in Cartesian coor- 
dinates. The main reason for using these coordinates is to allow posterior 
comparison of our work with other literature results. We demonstrate that 
considering a perfect fluid, there is no solution that incorporates rotation and 
expansion simultaneously that is in accord with Ozsvath [34]. By making the 
rotational parameters equals to zero, we obtain the Friedmann equations of 
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FLRW for a spatially flat universe. 

We also apply the Hamiltonian formulation and field equations of the 
TEGR implemented by Maluf [15], [16] to find the energy-momentum (grav- 
itational field plus matter) and gravitational angular momentum densities 
in the Godel-Obukhov universe, irrespective of the equations of the state of 
the cosmic fluid. Our result for total energy density is in accord with that 
obtained by Rybnickova, differing by a constant factor, where the Komar su- 
perpotential was used [11]. Our result for gravitational angular momentum 
disagree with that obtained by Dabrowski & Garecki using the gravitational 
angular momentum pseudotensor of Bergman- Thomson [10]. 

The article is organized as follows. In Section 2 we review the Lagrangian 
and Hamiltonian formulation of the TEGR. In Section 3, using the field 
equations of the TEGR, we find the teleparallel version of the Priedmann 
equations. In Section 4, we calculate the total energy of the Godel-Obukhov 
universe and compare it with those obtained from the pseudotensors. In 
Section 5, we find the total three-momentum of the universe. In Section 6, 
we also obtain the gravitational angular momentum density of this Godel- 
type universe. Finally, in Section 7, we present our conclusions. 

The notation is the following: space-time indices ... and global S0(3, 
1) indices a.b,... run from to 3. Time and space indices are indicated by 
H = 0,a = (0), (i). The fiat, Minkowski space-time metric tensor raises and 
lowers tetrad indices and is fixed by rjab = eaue^vQ^^ = ( — I- ++)• The 
determinant of the tetrad field is represented by e = det{e"' ^). We use units 
in which c = 1, where c is the light speed. 

2 The Hamiltonian constraints equations as 
an energy and gravitational angular mo- 
mentum equations 

We will briefiy recall both the Lagrangian and Hamiltonian formulations of 
the TEGR. The Lagrangian density for the gravitational field in the TEGR 
[35] with the cosmological constant A is given by 
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= -fc'eE"'"=r„bc - - 2eA;'A, (1) 

where k' — l/(167rG), G is the Newtonian gravitational constant, and stands 
for the Lagrangian density for the matter fields. As usual, tetrad fields 
convert space-time into Lorentz indices and vice-versa. The tensor E"**^ is 
defined by 

^ 2 

and — T"„^.Thc quadratic combination H'^'^Tabc is proportional to the 
scalar curvature i?(e), except for a total divergence. The field equations for 
the tetrad field read 



e„Ae6^9,(eE^^^) - e^E'''^ ^T^.^ - ^e^^T^edE''^'^) + ^ee„^A = ^ei;^ . (3) 

where eT^^ = bLul^e.""^- It is possible to prove by explicit calculations that 
the left hand side of Eq. (3) is exactly given by 

i e \RaM - ^e„/,i?(e) + e„^A| , (4) 

aand thus it follows that the field equations arising from the variation in L 
with respect to e"- ^ are strictly equivalent to Einstein's equations in tetrad 
form. 

The field equations (3) may be rewritten in the form 

(eE«^^) = ^ee% (t^'^ + f^'') , (5) 

where 

t^f^ = k' (4E^'=^r,/ - g^'^T.^'^'^n^d) , (6) 

and 

fx^ = T^^ + 2fc'^^^A, (7) 

are interpreted as the gravitational energy-momentum tensor [36] and the 
matter energy-momentum tensor, respectively. 
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The Hamiltonian formulation of the TEGR is obtained by first estab- 
hshing the phase space variables. The Lagrangian density does not contain 
the time derivative of the tetrad component . Therefore this quantity will 
arise as a Lagrange multiplier [37]. The momentum canonically conjugated 
to is given by 11"* = 5L/5eai. The Hamiltonian formulation is obtained by 
rewriting the Lagrangian density in the form L = pq — H , in terms of Cai, 
H"' and Lagrange multipliers. The Legendre transform can be successfully 
carried out, and the final form of the Hamiltonian density is expressed as [15] 

H = eaoC'' + a,,r' + Pkr\ (8) 

plus a surface term. Here a^fc and fS^ are Lagrange multipliers that (after 
solving the field equations) arc identified as =l/2(Tjofc + T^Qi) and j3k = 
Took ■ C"", r*^' and T'^ are first class constraints. The Poisson bracket between 
any two field quantities F and G is given by 

X 7 " ■^V5e,i(a;)5H»(a;) 6U-'{x) 6ea^{x) J ' ^ ' 

The constraint is written as = — c^jH^' + Zi", where h°' is an intricate 
expression of the field variables. The integral form of the constraint equation 
— motivates the definition of the energy- momentum four- vector [15] 

P^^ - j SxdiW' , (10) 

where V is an arbitrary volume of the three-dimensional space. In the con- 
figuration space we have 

Yiai = _4A;'eS"°^ . (11) 

The emergence of total divergences in the form of scalar or vector densities 
is possible in the framework of theories constructed out of the torsion tensor. 
Metric theories of gravity do not share this feature. 

By making A = in Eq. (5) and identifying H"* on the left hand side of 
the latter, the integral form of Eq. (5) can be written as 

P"= / £xee\[f^' + f'"') . (12) 

This equation suggests that is now understood as the total, gravitational, 
and matter fields (plus a cosmological constant fiuid) energy-momentum [36]. 
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The spatial components P*^*) form a total three-momentum, while a temporal 
component P*-"-* is the total energy (gravitational field plus matter) [1]. 

It is important to rewrite the Hamiltonian density H in the most simple 
form. It is possible to simplify the constraints which may be rewritten as 
a single constraint F"^. The Hamiltonian density (8) may be written in the 
equivalent form [16] 

H = e^oC"^ + l\ai,r^\ (13) 

where \ab = —\a are Lagrange multipliers that arc identified as Aj^ = 
and Aofc = — A^o = l^k- The constraints F"** = — F**" embody both constraints 
F*'^ and F*^ by means of relations 

= e„^e^,'=F«^ (14) 

and 

pfe = = F''^ = e„°e6*^F"^ (15) 

It reads 



p6a _ ^ab ^ (^aOb __ ^Wa^ _ ^^g^ 

Similarly to the definition of P", the integral form of the constraint equa- 
tion F"^ = motivates the new definition of the space-time angular momen- 
tum. The equation F"^ = implies 

= -Ak'e (e"°^ - Y.^") . (17) 
Therefore Maluf defines [16] 

L''^ ^ J^d^xe%e\M^''', (18) 

where 

M"'' = e^^e^M'^'^ = -M''". (19) 
as the four-angular momentum of the gravitational field. 
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The quantities and L are separately invariant under general coordi- 
nate transformations of the three-dimensional space and under time reparametriza- 
tions, which is an expected feature since these definitions arise in the Hamil- 
tonian formulation of the theory. Moreover these quantities transform co- 
variantly under global S0(3,l) transformations. 



3 Teleparallel version of Godel-Obukhov equa- 
tions 

The line element of the Godel-Obukhov or Godel-type cosmological model is 
given by the interval [30] 

ds'' = -df + 2a{t)V^e"'''dtdy + a{tf (dx'' + ke'^'^'^dy'' + dz'') . (20) 

where m, a and k arc the constant parameters and a{t) is the time-dependent 
cosmological scale factor. Here, we adopted the signature ( — A general 
analysis [30] , [38] of the kinematic properties of such a space-time shows that 
the closed time-like curves are absent in this manifold when A; > (cr > 
by definition). We obtain the usual Godel metric [26] by making a{t) — 1, 
cr = 1, m = 1 and k = —1/2 in the expression (20) . 
Using the relations 

and 



(21) 



a set of tetrads fields that satisfy the metric is given by [39] 

/ 1 -a(t)Vo'e'"^ \ 

a{t) 

a(i)VFTae"^^ 

V a{t) j 



(22) 



(23) 



This set of tetrads fields yields the velocity field given by e(o)'* (t, y, z) — 
(1, 0, 0, 0). According to the physical interpretation of Eq. (23), the latter is 
adapted to static observers in space-time [16]. 
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Now, with the help of the inverse metric tensor g^'^, we can write the 
inverse tetrads 

(24) 



as 



/ 1 

l/a{t) 

73^ 








V 











a(t)\/fc+c 





\ 





l/a{t) J 



(25) 



where the determinant of e"^ is 



e = a{t)^Vk + ae' 



(26) 



Before solving the field equations, it is necessary to consider the material 
content of the universe. We restrict our consideration here to the stress- 
energy-momentum tensor of a perfect fluid [40] given by 



{p + p) u^Uf, + pg^i^, 



(27) 



where p = p{x) is the matter energy density, p is the matter pressure and 
are velocity components. Here, we denote the velocity with respect to the 
commoving matter by 

which in covariant notation is written as -u^ = (^^a^f . Thus, we can write 
the nonzero components of the stress-energy-momentum tensor T^^^, as 



Too = {£ + p)uoUo + PQoo = £, (28) 

To2 = {s + p)uoU2+pgo2^T2o^-aV^e"''^e, (29) 

Til = (e + p)uiui + pgu ^ pa'^ , (30) 

T22 = ie+p)u2U2+pg22^a^e^"'''[ea+p{a + k)], (31) 

T33 = (£ + p)u3U3 + pg33 = pa^- (32) 

Now, it is convenient rewrite the field equations of the TEGR (3) as 



+ 
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where: 

+\ {ri^'e'^'e^^T^,, - r^^^e^'^e'^^T,^,) . (34) 
The nonzero components of the torsion tensor T^ixv are given by 



7^(1)01 = 


-7^(1)10 




(35) 


7(2)02 = 


— 2^(2)20 




(36) 


T(0)02 = 


"2^(0)20 




(37) 


^(0)12 = 


-T{Q)21 




(38) 


7'(2)12 = 


—T{2)21 


= may/k + ae"^^, 


(39) 



remembering that the torsion components are antisymmetrical under the 
exchange of the two last indexes. 

After tedious but straightforward calculations, we obtain the nonzero 
components of the tensor E"^'^ 



5](o)(i)(o) ^ 

J](3){3)(1) ^ 


_5](o)(o)(i) 


m 

~ 2^' 
m 

~ 2^' 




(40) 
(41) 


5](0)(1)(2) ^ 


_5](0)(2)(1) 


m 




(42) 


5](1)(2)(0) ^ 


_5](1)(0)(2) 


m j 




(43) 


^4^V 


k + a 


5](2)(i)(o) = 




m 

^ ~4^ 


]lk + a' 


(44) 






d 

a 




(45) 


5](2)(2)(0) ^ 


_j.(2)(0)(2) 


a 
a 




(46) 
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5](3)(3)(o) ^ _5](3)(o)(3) = (47) 

a 



5](o)(o)(2) = ^(o)(2)(o) ^ " / ^ (48) 

aVk + a' ^ ' 



^ = -^./^_, (49) 

aVk+a ^ ' 



5](3){3)(2) ^ X](3)(2)(3) ^ Q / ^ /^QN 

a]lk + a' ^ ^ 

Next, we proceed to obtain the components of the field equations in the 
model of the Godel-Obukhov universe. We mention here only the details 
used to write the components a = (0), /i = of field equations. We divide 
the field equations into five parts. Thus we rewrite the field equations as 



Parti = e(^a)xewd, {eec^e/E^'^) 



e(o)2e(o)(o)ao (ee(2)^e(o)°SW(^)W) + 
e(o)2e(o)(o)9i(ee(2)^e(i)^E(o)(2)(^)) 

aa + 

J- ^ V 



a a 



y/k + cr J y/k + a 



Part 2 = -ee/r]^o)d^'""^n,o 
ee(3)3s(3)(3)(o)T(3)3o 



a^aae^^ 



^k + a 

Parts = ^ee^^ee^e/TfeA.S^'^'^ 



iee„,(2e(o)°e(i)^r(i)oiE(^)(°)(^) + 2e(o)°e(2)^T(,)o2S(^)(°)(2) + 
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2e(3)^e(o)%)3oS(^)(^)(°)) 
\^kTa 2 8^/kTa 



Part A = -ee(o)oA = ^ A 



1 z.^ 1 0^ a^^e^^^sV^Ta 

Substituting Eqs. (23), (25), (26), (28), and (35)-(50) into the above equa- 
tion, and summing the parts, we obtain 

2(7 (aa — a^\ 3d^ ( k \ m^k 2 a „ ^ /.-.n 

+ 7^ - „2.,. , -uj'-A^87rGe, (51) 



+ (J y y \k + aj {k + a) 

where the dot indicates the temporal derivative and uj is the global angular 
velocity of the universe [31] given by 



m 
2a 



a 



(k + a) 



1/2 

(52) 



Repeating the procedure for the components a = (1), fi = 0; a = (1), 
fi = 1; a = (0), /i = 2; a = (2), /i = 0; and a — (3), /i = 3; we obtain, 
respectively, the independent equations 

a 

ma , = ~ 0, (53) 

2ad + a^ f k \ , a n ^ /r 

-u-A^ -SnGp, (54) 



\k + a 

-u^-A^ SnGe, (55) 



2a'd + f k \ rn^k ^ 



\k ^ o j a"^ {k + a 
{da — d^) 



„ = (56) 
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Equation (53) reveals the impossibility of the simultaneous existence of 
rotation and expansion using a perfect fluid matter source for the Godel- 
Obukhov metric. 

The field equations (53)-(57) represent the Godcl-Obukhov universe in 
the context of the TEGR. Fixing m = 0, cr = and A; = 1, they reduce to 
the Friedmann equations for the fiat universe 

— -A = SttCs, 

2aa + a^ . ^ . ^, 
^ A = -STrGp. (58) 

Moreover, Eqs. (53)-(57) are according to the works of Krcchct & Panov 
[41], and Korothii & Obukhov [42], using the Einstein equations and a perfect 
fluid 



4 Total energy of the Godel-Obukhov model 

Let us now calculate the total energy of the Godel-Obukhov universe using 
the equations shown in Section 2. By making A = 0, and substituting Eq. 
(5) in (12) and using 

^a\v ^ s^^^e/ce", (59) 

we have 

P" = Ak' I d^xdi (e S«^'=eb°ec') . (60) 

As already observed, the temporal component represents the system en- 
ergy. Therefore the energy will be given by 

P(°) = Ak' f d'xdi (eE(°)^^e6°ee') . (61) 

Such a quantity can be written as 

P(°) = 4k' I d'xd, (ee(o)°e(i)is(o)(°)(i) -t- ee^2)\i)'^^'^^'^^''^) • (62) 
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By substituting Eqs. (25), (26), (40), and (42) in (62), we obtain 

= -A;'am2^^±l / d'xe"'^. (63) 
Simplifying the above expression, it follows that the energy density is given 

by 

(o)__^2A^ 

" IQnGVkT^'' ■ ^^^^ 

We note that E^^'' is negative and clearly diverges when integrated in the 
entire space. 

Wc can now compare our results with those obtained by Rybnickova [11] 
and Dabrowski & Garecki [10]. 

Rybnickova obtained the total energy density, using the Komar su- 
perpotential, as 



to — 



Dabrowski and Garecki used the pseudotensor of Einstein for the station- 
ary metric of the acausal Godel model and arrived at the null result for the 
total energy density [10]. They found a negative total energy density for the 
causal Godel model. Sharif [43] calculated the total energy density associated 
with a space-time homogeneous Godel-type metric by using Einstein and Pa- 
papetrou energy- momentum complexes. Sharif's results are not in agreement 
with our result and he found that the two definitions of energy-momentum 
complexes do not provide the same result for this type of metric. 

Our result and that obtained by Rybnickova present the same dependence 
on a{t) and x, but differ by a constant factor. Moreover, we sec that E^^^ = 0, 
in the case k = —l/2,m = a = 1, corresponding to the energy density of the 
Godel universe (rotation only). Then, our result is compatible with those of 
the works of Rybnickova and Dabrowski & Garecki for the acausal model. 
They also found a nonzero total density energy for the causal Godel model. 

In particular, fixing m = a = (without rotation) and A; = 1 all these 
the results are identical to zero 

£(°)=u;° = 0, 

recovering the result for the fiat FLRW universe. 

These results for the fiat universe were also found by Rosen [3], Cooper- 
stock [4], Johri [5], and Garecki [6]. 
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5 The total momentum of the Godel-Obukhov 
model 



Let us now consider the calculation of the total three-momentum (matter 
plus gravitational field) of the FLRW Universe. As seen in Section 2, it is 
noted that the total three-momentum is given by space components a — {1}, 
{2} and {3} of Eq. (60) . 

In order to obtain the space component a = {1} of the total momentum, 
we can write the quantity P^^^ as 

P« = Ak' f d=^xai(ee(o)°e(i)iE(^)(°)(^) + ee(2)°e(i)iE(i)(2)(^)). (65) 

By substituting (25), (26), (45), and (49) in the previous equation, we 
obtain 

pW^^K^t j^d'xe^^ (66) 
Thus, the spatial momentum density p*^^) it is given by 

pW = ^^^^e™^ (67) 
AnG^/k -\- a 

For the absolutely analogous calculation, we have 

p« = !!^^e"-, (68) 

= 0. (69) 

We can observe that Eq. (67) is valid for any fluid matter. However, it 
shows that considering a perfect fluid, for a purely rotational or expansion 
universe, the spatial momentum density p^^^ is zero. Equation (68) shows 
that the spatial momentum density p^^^ is zero only for an expansion universe. 
We can note again that p^^^ and p^^^ diverge when integrated in the entire 
space. 
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6 The Gravitational angular momentum of 
the Godel-Obukhov model 



Let us verify the expression of gravitational angular momentum (18). By 
making use of (19) and (17), we can write (18) in the form 

By making use of S"*^^ = Cc^S"'^^ and using the determinant of tetrads, 
Eq.(26), and the components of the tensor S"^'^ , (40)-(50), we find that 
the components of the three-angular momentum densities and 
vanish. The unique nonzero gravitational three-angular momentum density 
component is given by 

£(1)(2) ^ _£(2)(1) ^ ^^y ^mx^ ^7^^ 

The above result demonstrates that there is only a direction to the gravita- 
tional angular momentum, as expected. The component >C(^)(^) found here 
diverges when integrated in the entire space. 

The other components are all null except for 

£{o)(i) = _£{i)(o) ^ moV^Ta^^^^ ^^2^ 

which, although it represents the component of the gravitational center of 
mass moment, does not possess physical meaning [44]. 

Our result Eq.(71) when fixing k — —1/2, m — a — 1 is not in agree- 
ment with those obtained by Dabrowski & Garecki which used the gravita- 
tional angular momentum pseudotensor of Bergmann-Thomson to calculate 
the gravitational angular momentum density [10]. Indeed as they used the 
pseudotensor to calculate the gravitational angular momentum, their results 
are not coordinate invariant. 



7 Conclusions 

In this work, we show explicitly the equivalence among equations obtained 
with the TEGR and those obtained by the GR for the Godel-Obukhov met- 
ric and also calculated the total energy-momentum and gravitational angular 
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momentum densities with the use of tensorial expressions of the TERG, ir- 
respective of the equations of state of the cosmic fluid metric. In the case in 
which the rotation parameters equals to zero, we recovered the flat FLRW 
universe. 

Our result for the energy density presents the same dependence on a{t) 
and X as that found by Rybnickova. Both results diverge when integrated in 
the entire space. Fixing k = —1/2, m = a = 1 (parameters of the stationary 
Godel metric), the energy density vanishes in accordance with the results 
obtained by Dabrowski & Garecki. 

By analyzing the equations for the total three-momentum (matter plus 
gravitational field) we conclude that all these components vanish simultane- 
ously in the universe without rotation. 

As we used tetrads fields adapted to static observers in space-time, we 
found a unique nonzero component of the gravitational angular momentum 
density with physical meaning, C'^^^^'^\ which refiects the preferred direction 
(z-axis) related to the cosmic rotation. 

We finally conclude from this work that the TEGR obtained equivalent 
results to the GR with the great advantage of addressing covariantly the def- 
initions of quantities as energy-momentum and angular momentum tensors 
of the gravitational field. 

In order to find a solution that simultaneously permits rotation and ex- 
pansion, we will consider others types of fiuids as the Chaplygin gas. Efforts 
in this respect will be carried out. 
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